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We investigate the cosmological moduli problem by studying a modulus decay in detail and find that 
the branching ratio of the gravitino production is genericaUy of 0(0.01 — 1), which causes another 
cosmological disaster. Consequently, the cosmological moduli problem cannot be solved simply 
by making the modulus mass heavier than 100 TeV. We also illustrate our results by explicitly 
calculating the branching ratio into the gravitinos in the mixed modulus-anomaly /KKLT- and 
racetrack-type models. 



The cosmological moduli problem T| is one of the most 
challenging puzzles in particle physics and cosmology. In 
this letter, we show that the problem is even more diffi- 
cult than usually thought. 

In supergravity/superstring theories, genericaUy there 
exist moduli fields which have flat potentials and obtain 
masses from supersymmetry (SUSY) breaking and non- 
perturbative effects. During an inflationary period, a 
modulus field X is likely to develop a large expectation 
value. After the end of the inflation, it starts a coher- 
ent oscillation and soon dominates the energy density of 
the universe. Due to the interaction suppressed by the 
Planck scale Mp = 2.4 x 10^^ GeV, the decay rate of the 
modulus X is extremely small: 
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which leads to an onset of a radiation-dominated uni- 
verse with a very low temperature: 



Tx = {n^g.m-^^WMpTx 

~ 5.5 X 10-^ MeV • f ^^"j 
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Here, c is an order one coefficient and is the effec- 
tive number of massless degrees of freedom. This is cos- 
mologically unacceptable because a successful big-bang 
nucleosynthesis (BBN) requires that the (last) radiation- 
dominated universe starts with temperature higher than 
~ 5 MeV 0. 

As is clear from Eq. (01, a simple solution would be to 
assume that the modulus X is ultra heavy ^: 



mx > 100 TcV 



Tx > O(MeV) 
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Actually, there have been proposed scenarios with such 
a large modulus mass (cf. [a, IE SB 01) ■ 

However, there exists yet another serious cosmological 
obstacle even for heavy moduli scenarios. The new prob- 
lem is caused by the gravitino which is produced by the 
modulus decay. Indeed, as we will show, the branching 



ratio of the modulus decay into the gravitino is generi- 
caUy quite large 
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Br(A: -> gravitino) = 0(0.01 - 1) 
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which causes serious problems after the modulus decay. 
We call this problem the moduli-induced gravitino prob- 
lem. 

The gravitino production via modulus decay and its 
cosmological implications have been previously discussed 
in Refs. (lOt [llf , taking Br{X gravitino) < 1. The 
main purpose of this letter is to show that Eq. holds 
in a generic setup, and to emphasize how disastrous its 
consequences are. We also exemplify explicit results in 
the mixed modulus-anomaly/KKLT mediation [EH ^^nd 
in the racetrack Q setups. 

Let us first estimate the branching ratio of a modu- 
lus decay into gravitino(s). We consider a heavy modu- 
lus scenario, mx ^ 100 TeV [cf. Eq. ©]. On the other 
hand, the gravitino is likely to be (much) lighter than 
100 TeV, because too large gravitino mass requires a fine- 
tuning in the Higgs sector due to the anomaly-mediated 
effects. Thus, we assume mx 3> W3/2 hereafter. After 
choosing the unitary gauge in the Einstein frame, where 
the goldstino component is completely absorbed into the 
gravitino under the super-Higgs mechanism |12}|. the rel- 
evant gravitino-modulus couplings are [TollllLll3l| 

e-'C = -^^'^P^ {GxdpX-GxdpX)^^,j,i^. 

-^e<=/2 {GxX + GxX) [7^ 7''] V., (5) 

where ip^ is the gravitino field. Here and hereafter, we 
set Mp = 1 unless explicitly written. The function G = 
if -I- In is the total Kahler potential, and Gi denotes 
a derivative of G with respect to the field i. The real and 
imaginary components of X ~ {Xp -f iXi)/\/2 have the 
decay rate ^, 
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See Refs. 



for other solutions. 



^ The branching ratio of the single gravitino production is 
suppressed due to the phase space; Br(X ^V'3/2) ^ 

(mg/a/mx)^. 
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in the limit of mx ^ 'ni?,/2 after the canonical nor- 
malization, X = y/Qxx-^^ where gij* = Kij- is the 
Kahler metric. We emphasize that the gravitino mass 
in the denominator arises from the longitudinal compo- 
nent of the gravitino, which corresponds to the goldstino, 
V'p ^ (Pa'/"^3/2)xgS) in the massless limit. This factor 
magnifies the decay rate by a factor of ~_(mx/?7T-3/2)^5 
compared to the previous results 0, 0, Q . 

The auxiliary field of the modulus, Gx, represents 
its fractional contribution to the total amount of the 
SUSY breaking. It may be small enough to suppress 
the gravitino-production rate. However, we can see that 
in the framework of the 4D N = 1 supergravity, Gx is 
generically bounded below, Gx ^ °'^'ij2l'^x = . The 
scalar potential is given by F = e*^(G*Gi — 3) and the 
(almost) vanishing vacuum energy demands that at least 
one of the auxiliary fields, Gi, should take G*Gi ^ 0(1), 
where G* — g^^ Gj* . If the modulus field plays the role, 
it dominantly decays into the gravitino [cf. Eqs.QJ and 
and hence clearly i?3/2 — 1. Instead, by introducing 
a hidden sector field Z with Gz ^ 0(1), the modulus- 
auxiliary field Gx can be small ^ . In the following dis- 
cussion, we assume that the modulus field does not decay 
into the hidden-sector field for simplicity. In order to see 
how small Gx can be, we investigate the condition to 
minimize the potential, Vx = 0, leading to 



G^^JxGx + G'-^JxGz + Gx = 0, 



(7) 



where ViG, 



G, 



rjj-Gfc with the connection, T^^ 



9^^ Uii'j- We assume that the Kahler potential includes 
no enhancement factor, and especially gxx^ 9zz ~ 
1 and Kx ^ 1- The first term in Eq. Q) then 
becomes ~ G^ R, because mx/m^/2 ~ |VxGx| + 
0{Gi). Barring cancellation, the second term is given by 
max [Wz/W, Wxz/W, Kxz, ^xz] ' which is 0(1) unless 
the hidden sector takes a special form. Thus, we arrive 
at \Gx\ - R'^ for \g^^\ < R'^ and \Gx\ > R'^ for 
^ i?"^. As a result, Eq. ® becomes 
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where we define Gx — k/R — Hi{m^/2/^x) with |k| ;C 1- 
Note that the above discussion is valid for any value of 
mx , as long as i? S> 1 is satisfied. 

The modulus field also decays into radiation, that is, 
the standard model (SM) particles and their superpart- 
ners. The relevant interactions of X then stem from the 
dilatonic couphng with the gauge sector, / (f6XWW, 



The potential may be up-lifted by adding the D-term, V^^K 
One can show that \Gx\ ^ also in this case, as far as 



V'_(,°)/y(D) 
Eq. Q. 



0(1), which corresponds to the second term in 



leading to 
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after the canonical normalization, where F^i, and A are 
the field strength of the gauge boson and its superpart- 
ner with a generator index a of the corresponding gauge 
symmetry. The chiral projection operators are defined 
as Vr^l = (1 ± 75)/2. We notice 6*^/^10^1 is approxi- 
mately given by the modulus mass for i? ^ 1. Then the 
decay rate is 
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for SU{3)c X SU{2)l x U{1)y, corresponding to c = 0(1). 
We notice that a half of the decay rate comes from the 
channel of the gaugino production, in contrast to the 
results of Refs. |lll IT3|. From Eqs. ijHl and l(T?Hl . we 
obtain the branching ratio of the gravitino production as 
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It is important to note that the production rate of the 
gravitino channel is one of the dominant processes in the 
modulus decay. Actually, it becomes an order of 0.01 — 1. 

The other decay processes are suppressed, except for 
a possible decay into the Iliggs(-ino) 14]. However, this 
decay channel is model-dependent, and does not change 
the above result c = 0(1). 

Let us show how such large branching ratio into the 
gravitinos jeopardizes the success of the standard big 
bang cosmology. We consider the constraints from (i) 
the speed-up effect, (ii) BBN, and (in) the lightest SUSY 
particle (LSP) abundance. First let us consider the so- 
called speed-up effect which modifies mostly the ^He 
abundance. This sets a bound on the abundance of the 
gravitinos with 7713/2 ^ 20 TeV (cf. 15]), since they de- 
cay after the neutron-proton transformation decouples. 
The observational data put an upper bound on the ratio 
of the energy densities of the gravitinos to the standard 
model particles at the BBN epoch as 
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In our heavy moduli scenario, however, the ratio is 
bounded below: /3/2 > B^^2/{^~ B3/2), where the equal- 
ity holds if the gravitinos are still relativistic at the BBN 
epoch and most of the superpartner of the SM particles 
directly produced from the modulus decay soon annihi- 
late into the SM particles. Thus we obtain 



B3/2 < 0.2, 



(13) 
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irrespective of whether the gravitino is stable or unstable. 
For B3/2 > 0.2, the gravitinos from the modulus decay al- 
ways upset the standard BBN, as long as m3/2 ^ 20 TeV. 

Now we discuss the cases of the stable and unstable 
gravitinos separately. First we take up the unstable grav- 
itinos, which is the case if the gravitino mass is heavier 
than the LSP mass, 7713/2 > trlsp- The gravitino-to- 
entropy ratio is given by 
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The BBN severely constrains ¥3/2 [isl Il7|. Even if we 
adopt the recent analysis on ^He abundance which has 
taken account of possible large systematic error [l8l |. 
13/2 must be smaller than 2 x 10"-^^ at 95% C.L. for 
m3/2 — 30 TeV 0, and the bound becomes much sev- 
erer for smaller 70,3/2. Therefore the branching ratio into 
the gravitinos must be extremely small: 
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FIG. 1: The cosmological bounds on 77x3/2 and -B3/2. Shaded 
regions are excluded by cosmological arguments. See the 
text for details. The horizontal dashed line denotes the 
BBN bound from the stau NLSP decay into gravitinos for 
777NLSP ~ 100 GeV. We have chosen nix ~ TeV and 
c = 1 as representative values. The bounds become severer 
for larger mx- 
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for 777LSP < 7773/2 < 30 TeV. Here e < 1 parameterizes 
the BBN bound: e = 1 for 7773/2 ~ 30 TeV, and 10"^ < 
e <C 1 for 7713/2 < 30 TeV. In addition, the abundance of 
the LSPs from the gravitino decay is (cf. (l^) 



Ilsp 



3/27 



45 



1 



StT^^, MpT'i/2 {<Ja 



where {uannv) is the thermally averaged annihilation 
cross section of the LSP, and T3/2 is the decay temper- 
ature of the gravitino. Since the LSP abundance must 
be smaller than the dark matter abundance, we have an- 
other constraint on -B3/2: 



(^100GeV\ 
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(17) 

for 7773/2 > TTiLSP- Here JIdm is the density parame- 
ter of the dark matter, h is the present Hubble param- 
eter in units of lOOkm/sec/Mpc, and we have assumed 
{(Tannv) < 10"^ GcV^^ and 7773/2 < 100 TeV. We can see 
that (|13|l . H15|) and (|17|l rule out the unstable gravitinos, 
unless i33/2 is extraordinarily small in spite of our result 
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Next we consider the stable gravitinos, which is the 
case if the gravitino is the LSP. A constraint then 
comes from the requirement that the gravitino abun- 
dance should not exceed the dark matter abundance, and 
we only have to replace ttilsp with 7773/2 in H17(l : 
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for the stable gravitino. Furthermore, the bound on the 
gravitino abundance is severer by an order of magni- 
tude for 7773/2 < 100 MeV(777x/103 TeV)-l/^ due to 
the present large free-streaming velocity [23 . See Fig. ^ 
From (|13|l and l|18|) . we conclude that large B3/2 encoun- 
ters trouble even for the stable gravitinos. 

Lastly we comment on the lightest superpartner of the 
SM particles, denoted by x, produced from the mod- 
ulus decay through the gauginos. From the discussion 
above, the number of x generated from the decay of 
unit quantum of the modulus field is order unity. If x 
is the LSP and stable, it must be electrically neutral. To 
satisfy J^lsp^^ < 0.13. the pair annihilation cross sec- 
tion must be large |l4lll9| . which in turn constrains the 
mass spectrum of the SUSY particles. In the case of the 
gravitino LSP, there is a strict BBN bound on the abun- 
dance and lifetime of the next-to-LSP (NLSP) %. For 
a stau NLSP f , we find upper bounds on the gravitino 
- 1) GeV for mx = 10^ TeV and 
TeV For a neutralino NLSP, the 



mass, 



< (0.3 

777 f 5: 1 



^773/2 
100 GeV < _ 
bound becomes severer. 

We summarize the bounds considered above in Fig. ^ 
from which we can see how serious the cosmological mod- 
uli problem becomes as a result of our finding that B3/2 
should be 0(0.01 — 1). Note that the bounds on ^3/2 [cf. 
(|15|l . (|17|l and (|18|l ] become severer for heavier modulus 
mass and larger c, although the dependence is weak. 



dWe have used ~ (45/87r2g.)i/2 (MpTy (aa„„-!;»-l [T| 

with ((Tannv) < IQ-^mZ^ BBN bounds from 

I15I with t's hadronic branching in |22| . and Tf ~ 
6sec(mf/100 GeV)-'^ (7713/2/10 MeV)^ . The bound from ^3/2 
(from -f ) < Qdm is slightly weaker than the BBN bound. 
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Let us explicitly calculate the branching fraction of 
the gravitino in some models. The relevant terms in the 
Lagrangian of the mixed modulus-anomaly/KKLT me- 
diation model 0, Q are given by 

/ = -3{X + X)'"/^ + Cz{X + xy\Z\^, (19) 



W ^ Wq + ae 
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where the Kahler potential is defined as / = — 3e~^/^. 
Here coefficients b and Cz are real with hXji 3> 1, while 
a and Wq are complex. Then the modulus mass is 
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(21) 



The gravitino mass TO3/2 is much larger than the soft 
mass scale |2j- On the other hand, Gx is 
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where we assumed (Z) <C 1 and the vanishing cosmolog- 
ical constant. Thus the branching ratio becomes 



BiiXi ^ 2^3/2) 



(2n + 30V216 
l + (2n + 30V216' 



(23) 



In the original KKLT model [g with the lift-up potential 
due to D3 brane, which is realized by choosing n = 3 and 
I = 0, the branching ratio is -B3/2 — 0.14. 

In the racetrack-type models ^] , the Kahler potential 
is the same as Eq. H19|l but the superpotential is 

W = aie-''!^ + 026-''^^ + VFHiddcn(^), (24) 

where a^'s are complex and 6i's are real numbers. Here 
biXfi 1. The modulus mass and Gx are 
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Therefore the branching ratio becomes the same as 
Eq. For n = 1 and ^ = 0, one obtains B3/2 ~ 0.018. 

In summary, we have shown that the branching ra- 
tio of the modulus decay into gravitinos is 0(0.01 — 1), 
and hence heavy modulus scenario is plagued with the 
moduli-induced gravitino problem. Here let us comment 
on possible ways out. Unless either i33/2 or the grav- 
itino mass is extremely small, there is no way to cir- 
cumvent the problem other than introducing something 
that dilutes the modulus field and the subsequently pro- 
duced gravitinos and (N)LSPs. One of the candidates 
is the thermal inflation 23]. Another is Q-balls in 
the Affieck-Dine mechanism |23. In the latter case, the 
Q-balls can not only dilute the unwanted gravitinos but 
also generate the baryon asymmetry successfully. De- 
tailed discussion is beyond the scope of this letter and 
will be presented elsewhere p^ . 

We thank Wilfried Buchmiiller for comments. M.E. is 
grateful to Shuntaro Nakamura and Masahiro Yamaguchi 
for private communication. M.E. and F.T. would like to 
thank the Japan Society for Promotion of Science for 
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Note Added: Shortly after our letter, Ref. |23| ap- 
peared, pointing out the same problem. They also dis- 
cussed the gravitino much heavier than lOOTeV as a pos- 
sible solution to the moduli-induced gravitino problem. 
More recently, Ref. has pointed out Gx can be much 
smaller than in some situations. However, for gen- 
eral non-rcnormalizablc couplings in K , our estimate on 
the gravitino production still holds. See [l^ for details. 
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